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ABSTRACT 

The redshift-space distortion (RSD) of galaxies surrounding massive clusters is 
emerging as a promising testbed for theories of modified gravity. Conventional 
applications of this method rely upon the assumption that the velocity field in the 
cluster environment is uniquely determined by the cluster mass profile. Yet, real 
dark matter halos in N-body simulations are known to violate the assumption 
that virial mass determines the configuration space distribution, an effect known 
as assembly bias. In this Letter, I show that assembly bias in simulated dark 
matter halos also manifests in velocity space. In the 1 — 10 Mpc environment 
surrounding a cluster, high-concentration “tracer” halos exhibit a 10 — 20% larger 
pairwise-velocity dispersion profile relative to low-concentration tracer halos of 
the same mass. This difference is comparable to the size of the RSD signal 
predicted by f(R) models designed to account for the cosmic acceleration. I use 
the age matching technique to study how color-selection effects may influence the 
cluster RSD signal, finding a ~ 10% effect due to redder satellites preferentially 
occupying higher mass halos, and a ~ 5% effect due to assembly-biased colors 
of centrals. In order to use cluster RSD measurements to robustly constrain 
modified gravity, we likely will need to develop empirical galaxy formation models 
more sophisticated than any in the current literature. 

Key words: cosmology: theory — dark matter — galaxies: halos — galaxies: 
evolution - large-scale structure of universe 


1 INTRODUCTION 


The observed acceleration of the universe at late cosmo¬ 
logical times has been one of the major puzzles of fun¬ 
damental physics since its discovery. In the standard cos¬ 
mological model, ACDM, the cosmic acceleration is ac¬ 
counted for by dark energy, a mysterious, perfectly ho¬ 
mogeneous cosmic fluid with a constant equation of state 
parameter w = — 1. Alternatively, the cosmic acceleration 
could instead be driven by modifications to general rel¬ 
ativity on cosmological scales. The quest to uncover the 
physical nature of the acceleration is a primary science 
target of numerous current and planned galaxy surveys. 


Observations of the redshift-space distortions of 


els of modified gravity (IZhang et al. 2007: Linder 

2008; 

Reves et al. 

2010i: Dennings et al.l 20121: Reid et al 

2014; 

Beutler et al. 2014!) . In narticular. the velocity structure 


surrounding massive galaxy clusters is emerging as one of 
the most sensitive large-scale structure prob e s of depar¬ 
tures from general relativit y (ISchmidtl |20l d: iLam et al.l 
120121 . l2013l : IZu fe WeinbeiJl2013l : Zu et al. 2014 ). Under¬ 


standing the infall region of a massive cluster requires 
making predictions for the velocity field deep into the 
nonlinear regime, for which th e halo model i s the premier 


theoretical tool (iSeliak 

l200ll; IWhitd l200ll; ISheth et al. 

l200ll: iKang et al.1 20oJ 

Tinker et al. 20061:iTinker 2007). 


The baseline proving grounds for any precision model 
of velocity structure are cosmological N-body simula¬ 
tions. In particular, for models attempting to predict 
the velocity field on the highly nonlinear scales of ~ 
1 — lOMpc, it is commonplace to evaluate a model’s suc¬ 
cess by its ability to accurately predict the RSD signal 
of interest exhibited by host halos selected by mass. For 
example, a velocity structure model may be deemed ac¬ 
curate if the distribution of pairwise-velocities of dark 
matter halos is accurately predicted as a function the 
masses and spatial separation of the halo pairs. 

There are many reasons why the above class of cal¬ 
ibration tests can only be considered a prerequisite to 
the development of a complete model of RSD suitable 
to potentially falsify general relativity with galaxy red- 
shift surveys. Galaxies are biased tracers of dark matter 
halos, and the host halo pairwise-velocity field depends 
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on halo mass (lTinkerll2007l 'l. Another contribution to the 
pairwise-velocity field comes from the velocity disper¬ 
sion of satellite gal axies orbiting within the po tential well 
of their host halo (IHikage fc Yamamotcll2013h. an effect 
with a pronounced halo mass -dependence dMore et al.l 
120091 : IWoitak fe Mamonl l2013h . These effects naturally 
vary with the selection function of the galaxy sample, 
and so any model for the velocity field must be flexible 
enough to encapsulate variations due to the manner by 
which galaxies weight halos by mass. 

Although mass is the dominant variable determining 
the phase space distribution of dark matter halos, it is by 
now firmly established that knowledge of halo mass alone 
is insufficient to characterize halo clustering statistics, an 
effect that goes by the catch-all term halo assembly bias 
|Gao et al.||2005|: Wechsler_e£_alJ 200d; Wang_et_ay 200ft 


Gao & White 2007: Dalai et al. 2008: Lacerna & Padilla 


201ll l. In addition to configuration space-based statistics 


such as two-point clustering, the i ntra-halo velocity struc¬ 
ture o f both dark matter partinl«G[Fjdtenbache7|&3Vhi2 
l20ld l and orbiting subhalos dFaltenbacherl2Qlol 'l has been 
shown to depend upon host halo properties besides mass. 
However, halo assembly bias in the two-halo regime of the 
velocity field has received considerably less attention. 

If the statistical connection between halos and their 
resident galaxies depends upon properties besides halo 
mass, then the galaxy population will inherit the assem¬ 
bly bias of its parent halos. Such a population is said to 
ex hibit galaxy as s embly bias. It has recently been shown 
in IZentner et alj (12013 ) that galaxy assembly bias can 
lead to systematic errors that seriously threaten the pro¬ 
gram to use conventional implementations of the halo 
model to interpret obser vations of projected two-point 
clustering. The findings in lZentner et al.l (I2013l f are based 
on an empirical model of galaxy formation (described in 
m, though similar conclusi ons apply to results based on 
semi- analytic models fe.g.. ICroton et al] 120071 : IZu et all 

l2008h . 


Motivated by these findings, in this Letter I take a 
first look at how color-dependence of the galaxy selec¬ 
tion function may influence the redshift-space distortion 
signal observed in the neighborhood of massive galaxy 
clusters. I describe the simulations and mock catalogs I 
use throughout the paper in § [5] I present my results in 
§[3] and discuss their implications in §[4] 


2 SIMULATION & MOCK CATALOGS 

To study the velocity field surrounding massive ha- 
los, I use the co llisionless N-body Multidark simulation 
llRiebe et al.ll201lh . The Mult idark cosmological pa rame- 
ters are based on WMAP5 dKomatsu et al.ll2009ll . and 
the simulation was run with 2048 3 particles of mass 
m p « 8.7 x 10 9 Mq//i in a lGpc//i periodic box. In all the 
results in this paper, I use ROCK STAR-based halo cat¬ 
alogs llBehroozi et al.l 1201 ll . l2013l l , publicly available at 
http://hipacc.uese.edu/Bolshoi/MergerTrees.html 
In addition to the Multidark halo catalogs, I also 
study cluster velocity structure usi ng a mock galaxy cat¬ 
alog based on abundance matching (iKravtsov et ai1 l2004 ; 
IConrov et al.ll2006l : iBehroozi et al.ll2010h and age match¬ 


ing dHearin fe Watsonll2013l 'l. In particular, I use the cat¬ 
alog based on M r luminosity and g — r color, publicly 
available at http://logrus.uchicago.edu/~aphearin The 
value of using the age matching catalogs is that they pro¬ 
vide a reasonably realistic represe ntation of the galaxy 

12 : 


2012 


distribution in the local universe dRedd ick et al. 
iHearin et al.l 120121 . 120131 : I Watson et aid 1201411. and they 
exhib it strong signatures of assembly bias ( Zentner et al.l 
120131 '). In order to parse effects on the velocity field that 
are due to assembly bias from effects that are purely due 
to the manner by which halos are weighted by galaxies 
as a function of halo mass, I additionally use a mock 
galaxy catalog in which M v n is the only halo property 
governing the colors of the halo’s galaxies. I construct 
this “no-assembly-bias” mock catalog simply by shuffling 
the colors between mock galaxies occupying halos of a 
similar mass, using mass bins O.ldex in width, and sepa¬ 
rately shuffling the colors of centrals and satellites. 


3 RESULTS 

3.1 Overview of Methodology 

In both sections below, I study the pairwise-velocity dis¬ 
persion of dark matter halos, denoted by a(R\Mi, M2), 
where Mi and M2 are the masses of the halo pairs and R 
is their spatial separation. Motivated b y the program t o 
constrain modified gravity described in I Schmidt (l201(i) . 
I will focus on the subspace of information provided by 
the special case where 2 Mpc /h < R < 10 Mpc/h and 
one of the halo pairs in a(R\M\, M2) pertains to a cluster 
whose mass M 2 ~ 10 1A Mq/ h is known through indepen¬ 
dent means, for example through gravitational lensing. 

In all that follows, I keep the cluster sample fixed and 
study how the cluster velocity dispersion profile depends 
on the properties of the “tracer” Mi halos. In H3.2I I show 
that not only does velocity structure depend on tracer 
mass Mi, but has an additional dependence on tracer 
concentration: 

<j(R|Mi;M 2 = 10 14 M o ) ^a(R|Mi,ci;M 2 = 10 14 M Q ), 

the manifestation of halo assembly bias in velocity space. 
Of course, we do not have direct observational access to 
Mi, and any tracer galaxy sample will occupy a range of 
halo masses. So in <13.31 1 use the age matching mocks to 
show how the galaxy occupation statistics of the tracer 
halos influence mock observations of the velocity disper¬ 
sion profile, again studying both Mi—dependent effects, 
as well as effects from assembly bias. 


3.2 Dependence of the Cluster Velocity Field 
on the Tracer Halo Population 

In this section, I demonstrate that the velocity field sur¬ 
rounding massive cluster halos is sensitive to properties 
of the tracer halo population besides M v i r . I calculate 
the velocity structure surrounding Multidark host ha¬ 
los as follows. First, I randomly select a sample of 10 4 
cluster-sized host halos with mass 10 13,9 Mq/L < M v i r < 
10 141 Mq//i, which has a median R v n « 950kpc /h. I se¬ 
lect a second sample of host halos of mass 10 12 , 5 Mq//i < 
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Figure 1 . Pairwise radial velocity dispersion profile of 
tracer halos surrounding cluster-sized halos, plotted as a func¬ 
tion of cluster-centric distance. Differences between the solid 
and dashed black curves show the dependence of the signal 
on the mass of the tracer population. Red and blue curves 
show results for a tracer sub-population with high- and low- 
concentrations for their mass, respectively. The difference be¬ 
tween the red and blue curves demonstrates the manifestation 
of assembly bias in velocity space. Differences in are at the 
~ 10% level, of order the difference generated by departures 
from general relativity predicted by /( R.) models designed to 
account for the cosmic acceleration. 


M v i r < 10 13 ' 5 Mq//i, which I use as tracers of the veloc¬ 
ity field surrounding the cluster-sized halos. In spherical 
shells surrounding each cluster halo, I calculate the radial 
component of the relative peculiar velocity between the 
cluster halo and every tracer halo in the shell, denoted by 
I4 a d- I then stack all the values of Kad for every tracer of 
every cluster, compute the dispersion er(? d , and plot the 
result with the black curve in Figure [T] as a function of 
the radial distance to the mid-point of the spherical shell. 

With the red and blue curves in Figure [T] I ad¬ 
ditionally show the velocity structure traced by high- 
and low-concentration tracer halos in the mass range 
10 12 ’ 5 Mq//i < M v i r < 10 13 Mg/ft. To obtain these 
concentration-selected tracer sub-populations, I bin the 
tracer halos by mass using 0.02dex-width bins, and split 
each bin by its median value of NFW concentration, so 
that the two sub-populations have the same mass func¬ 
tion. 

For R > 2 Mpc /h, radial velocity dispersion a™ d 
decreases with decreasing cluster-centric distance: the 
relative velocity of halos becomes more coherent as 
the tracers share an increasing number of larg e -scale 
density modes sou rcing their velocities dFisherl Il995l : 
iReid &; Whitdl201ll '). The dependence of ay d on tracer- 
mass reflects the fact that higher mass tracers are less 
susceptible to tidal forces, so higher mass tracers flow 
more coherently with the large-scale velocity field (see 



Rp[Mpc/h] 


Figure 2. Pairwise line-of-sight velocity dispersion profile of 
mock galaxies surrounding clusters, plotted as a function of 
the projected cluster-centric distance. Lower, thin curves ex¬ 
clude satellites from the tracer population, upper, thick curves 
include the satellite contribution. Black curves pertain to a 
tracer population of all galaxies; red curves show results af¬ 
ter applying an additional red sequence cut on the tracers. 
Solid curves use the age matching mock, dashed curves its 
counterpart in which assembly bias effects have been erased. 
Differences between the red and black corresponding curves 
demonstrate the significance of color-selection effects on the 
cluster RSD signal. Differences between the top and bottom 
trio of curves indicate the critical role of accurate modeling 
of satellite velocities. Differences between the red dashed and 
red solid curves show the influence of assembly bias. 


also the top panels of Figure 6 of lTinkeill2007l . for an al¬ 
ternative demonstration of this point). The most striking 
effect shown in Figure |T| is the trend in tracer concen¬ 
tration. Halos with above-average concentration for their 
mass exhibit ~ 10 — 20% higher <j T y d relative to below- 
average concentration tracers. 

The trend in cry ' d with tracer concentration can be 
understood in terms of the well-established trend of two- 
point spatial clustering with concentration. Higher con¬ 
centration halos cluster more strongly relative to lower 
concentration halos of the same massQ Thus at fixed 
mass, high-concentration halos a re found in pr eferentially 
over-dense regions. As shown in lTinkeil (|2007h , the halo- 
halo pairwise-velocity distribution depends upon large- 
scale density, even when the masses of the halo pairs are 
held fixed. Thus sub-selecting tracer halos with high (low) 
concentration effectively selects pairs in high (low) den- 


1 Note that this statement only applies to the mass range rel¬ 
evant to the tracer population shown here; this trend reverses 
for h alos with masses M v j r M*, the halo model collapse 
mass (IDalai et al.1120081 ). 
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sity environments, naturally giving rise to the character 
of the trend seen in Figure [T| 


3.3 Dependence of the Cluster Velocity Field 
on the Tracer Galaxy Population 

In this section, I use the age matching mock galaxy cat¬ 
alog to translate the results of § 13.21 into an observa¬ 
tional context. For definiteness, I wi ll focus exc l usivel y 
on the cluster RSD signal studied in lLam et all J2012I '). 
the line-of-sight velocity dispersion profile in the cluster 
environment. Note, though, that the galaxy occupation 
principles discussed below apply equally well to alterna¬ 
tive statistical quantifications of RSD. 

I begin by selecting all ~ 10 14 Mq//i clusters, de¬ 
fined as galaxy systems residing in host halos with mass 
10 13 ' 15 M@/h < M v i r < 10 14 ’ 25 Mq//i. Treating the sim¬ 
ulation z-axis as the line-of-sight, I place mock galaxies 
into redshift-space. For a given tracer population (de¬ 
scribed below), I select the tracers residing in a cylindri¬ 
cal annulus 4000 km/s in length, and compute the line- 
of-sight velocity of the tracers in the annulus relative to 
their associated cluster. Figure [2] shows a(? s , the relative 
line-of-sight velocity dispersion, plotted as a function of 
the distance to the mid-point of the cylindrical annulus. 

Figure [2] shows results for the above calculation per¬ 
taining to several different choices for tracer population 
of mock galaxies. In the bottom trio of thin curves, I 
consider an idealized observation in which central galaxy 
identification is perfect, and all satellites have been ex¬ 
cluded from the tracer population. In the top trio of thick 
curves, I consider the opposite extreme case where no at¬ 
tempt whatsoever is made to exclude the satellite con¬ 
tribution. In practice, if some BCG-selection and/or iso¬ 
lation criteria is applied to the tracer population, the 
measured values will lie between the top and bottom trio 
of curves. Black curves show results for the full popula¬ 
tion of M r < —19 tracers, while red curves pertain to 
tracers passing an additional red sequence cut, defined 
by {g — r) = 0.21 — 0.03M r . Solid red curves use colors 
taken directly from the age matching mock, while dashed 
red curves pertain to the second mock galaxy catalog in 
which assembly bias has been erased by the shuffling pro¬ 
cedure described in § [2] 

The first general trend to note is the change in the 
slope of all curves in Figure [2] relative to those in Fig¬ 
ure [0 This difference can be understood through the 
characteristic shape of the pai rwise-velocity dispe rsion 
profile as a function of scale IITink cr et al. 120061 . Fig¬ 
ure 4). On sufficiently small scales, the dispersion pro¬ 
file turns over and increases with decreasing pair sepa¬ 
ration due to the mutual gravitational influence of the 
halo pairs. The radial dispersion does not exhibit this 
turnover until R < 2Mpc//i, but line-of-sight smearing 
together with the Az cut result in the turnover occuring 
at R p « 10Mpc//i. 

Let us now unpack how the galaxy sample selection 
function can influence the cluster RSD signal by com¬ 
paring various pairs of curves in Figure [2] The simplest 
comparison to understand is the difference between the 
top and bottom trio of curves. In all cases, the intra-halo 


velocity dispersion of satellites significantly enhances cr(/ s 
on all scales, as should be expected. 

Next consider the signal exhibited by the idealized 
centrals-only samples shown with the lower trio of thin 
curves. At the fainter end of the luminosity function, age 
matching predicts that red centrals reside in slightly lower 
mass halos relative to blue centrals of the same bright¬ 
ness. Recall from Figure[T]that ay decreases with increas¬ 
ing halo mass. In the model with no assembly bias, this is 
the only operative trend for centrals, and so red centrals 
living in less massive halos exhibit a larger cry relative 
to all centrals: the thin, red dashed curve lies above the 
thin, solid black curve. In the age matching model that 
includes assembly bias, red centrals not only live in lower 
mass halos, but red centrals additionally occupy halos 
of above-average concentration. This assembly bias effect 
further boosts the velocity dispersion profile of red se¬ 
quence centrals: the thin, solid red curve shows a ~ 5% 
boost relative to the thin, dashed red curve. This assem¬ 
bly bias boost is simply the concentration trend shown 
with the red and blue curves Figure [I] as manifested in 
the age matching mock. 

Finally, consider the top trio of thick curves illus¬ 
trating results for tracer galaxies that include a satellite 
contribution, beginning by comparing the thick, dashed 
red and thick, solid black curves. This comparison demon¬ 
strates the effect that a red sequence cut has on ay for a 
galaxy population in which color is governed exclusively 
by virial mass. Any red sequence cut has two, related 
effects on the satellite population: 

(i) At fixed luminosity, redder populations of galaxies 
have a larger satellite fraction; 

(ii) redder populations of satellites are found in pref¬ 
erentially more massive host halos. 

Both of these “M v i r — only” selection effects work together 
so that a red sequence cut generically boosts ay s due 
to an increased contribution from the internal motion of 
satellites within their host halos. 

Finally, let us consider the role played by assembly 
bias in the tracer samples that include satellites by ex¬ 
amining the thick, solid red curve. In the galaxy sample 
illustrated by this curve, the color-selection effects on the 
satellite population discussed in the previous paragraph 
still apply, but there is an additional effect due to as¬ 
sembly bias in both the centrals and satellites. From the 
bottom trio of curves, we know that assembly bias in the 
central population will generically lead to a boost cr(? s af- 
ter applying the red sequence cut. However, as shown in 
IZentner et al.l d2013l f . age matching predicts much weaker 
assembly bias effects for satellites relative to centrals (see 
also McEwan et al. 2015, in prep). Thus although as¬ 
sembly bias is still operative in the sample that includes 
satellites, the “M v i r — only” effects dominate. For mea¬ 
surements made on tracer samples with no satellite ex¬ 
clusion criteria, age matching predicts the role played by 
assembly bias in color-selection effects on to be lim¬ 
ited to 1 — 2%, as shown by the difference between the 
thick, solid red and thick, solid dashed curves. 
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4 DISCUSSION 


There is rich information about our cosmology encoded 
in P(AV i 2 |-R, Mi, M 2 ), the pairwise-ve locity distribution 
of halos llTinkerll2007l : iReid et ~ai]l2014T ) . Recent work has 
focused on the special case where one of the halo pairs is a 
clust er with mass that is assumed to b e known by lensing 
(e.g., lLam et al ■ 2013; Zu et alj |20lj). Preliminary indi¬ 
cations are quite promising, as these authors have shown 
that cluster velocity structure can be modeled with rea¬ 
sonably good accuracy, and with a signal that varies at 
the 10 — 20%-level due to modified gravity effects. 

In practice, we must rely upon galaxy redshift sur¬ 
veys to provide our information about the velocity field, 
and so regardless of the approach to the RSD model¬ 
ing, the method must be able to account for the biased 
manner by which the tracer galaxies occupy halos. The 
results in A3.3I based on the age matching mocks show 
that the effects of the galaxy selection function on the 
cluster velocity field are roughy as large as the modified 
gravity signal itself. These findings imply that the pro¬ 
gram to place precision constraints on general relativity 
with cluster RSD will need to proceed hand-in-hand with 
realistic models connecting galaxies to dark matter halos. 

The success of age matching implies it is at least 
plausible that galaxy color correlates with halo assem¬ 
bly at fixed halo mass. However, traditional implemen¬ 
tations of empirical galaxy-h alo models such as the Halo 
Occupation Distribution (IBerlind fe WeinbereH2002r ) and 


B 

fj 2003) 


Conditional Luminosity Function (lYang et al.l 20031) ne¬ 
glect to account for assembly bias; such implementa¬ 
tions are unable even in principle to encapsulate the 
complexity shown in Fig. [2] The semi-ana lytic approach 
to the galaxy-halo co nnection (see, e.g., iBenso 3 120121 ; 
iHenriaues et al.|l2014l . and references therein) is not lim¬ 
ited by such simplifying assumptions, but at present this 
class of models is too computationally demanding to use 
in cosmological parameter likelihood analyses. In order to 
robustly exploit the cluster RSD signal that will be mea¬ 
sured with Stage IV dark energy experiments, we likely 
need to develop a new generation of galaxy-halo models 
more sophisticated than any in the current literature. 
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